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INTRODUCTION 

In  Ref.  1,  a preliminary  investigation  has  been  made  for  the 
reliability  prediction  of  composites  under  periodic  proof  tests 
in  service,  along  with  the  determination  of  the  optimal  periodic 
proof  test.  The  theory  of  periodic  proof  tests  and  the  reliability 
prediction  for  composites  has  been  based  on  a particular  residual 
strength  degradation  model  (Refs.  1-4).  It  is  the  purpose  of  this 
paper  to  verify  experimentally  the  theory  of  proof  tests  established 
in  Ref.  1 and  the  validity  of  the  residual  strength  degradation 
model  used. 

Because  of  the  importance  of  the  residual  strength  degradation 
model  to  the  reliability  prediction  of  composites,  a new  theoreti- 
cal derivation  of  such  a model,  based  on  the  assumption  that  the 
residual  strength  degrades  monotonically , is  presented.  The  experi- 
mental verification  of  such  a model  using  constant  amplitude  fatigue 
data  is  then  followed  in  Section  II. 

After  the  experimental  verification  of  the  residual  strength 
degradation  model,  two  series  of  experimental  data  under  periodic 
proof  tests  are  presented  in  Section  III.  The  theoretical  predic- 
tion of  the  residual  strength  degradation  under  both  cyclic  fatigue 
loading  and  periodic  proof  tests  is  derived  and  correlated  with 
experimental  data. 

It  is  shown  that  (i)  the  theoretical  residual  strength  degra- 
dation model  correlates  reasonably  well  with  the  experimental  data 
for  Gr/E  unnotched  laminates  under  constant  amplitude  fatigue  loading, 
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(ii)  the  theory  of  periodic  proof  tests  established  in  Ref.  1 
correlates  extremely  well  with  the  experimental  data,  (iii)  damage 
to  Gr/E  unnotched  laminates,  in  term  of  static  strength  and  fatigue 
strength,  due  to  repetitive  proof  tests  appears  to  be  negligible 
for  the  proof  load  level  up  to  901  of  the  mean  ultimate  strength, 
and  (iv)  the  residual  strength  degradation  model,  the  theory  of 
periodic  proof  tests,  and  the  methodology  for  reliability  predic- 
tion of  composites  under  periodic  proof  tests  in  service  esta- 
blished in  Ref.  1 can  be  used  with  confidence  for  Gr/E  unnotched 
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SECTION  II 


RESIDUAL  STRENGTH  DEGRADATION  MODEL 


(1)  Theoretical  Derivation 

Experimental  fatigue  data  (Refs.  2-3)  on  unnotched  composite 
laminates  indicate  that  the  residual  strength  R(n)  after  n fatigue 
cycles  is  a monotonical ly  decreasing  function  of  n.  Assuming  that 
the  slope  of  the  residual  strength  R(n)  is  inversely  proportional 
to  some  power  c-1  of  the  residual  strength  R(n)  itself,  one  obtains 
a rate  equation, 


dR(n)/dn  = -f  (S,u),e)/cR‘^‘ Mn) 


in  which  f(S,d),e)  is  a function  of  the  maximum  cyclic  stress  S, 
the  loading  frequency  w,  and  the  stress  ratio  e.  In  Eq.  1,  c is 
a constant.  For  the  sake  of  simplicity,  the  loading  frequency  u 
and  the  stress  ratio  e will  be  fixed  so  that 


f(S,o),e)  = f(S) 


Integration  of  Eq.  1 from  nQ  to  n^^  cycles  leads  to  the 
following  result. 


R^'Cnj)  = R^'Cnp)  - f(S)(nj-np) 


in  which  f(s)  is  a function  of  maximum  stress  S. 


With  ng  = 0,  nj^  = n,  Eq.  3 reduces  to 
R'^(n)  = r'^(O)  - f(S)n 


Eq.  4 expresses  the  residual  strength  at  n cycles,  R(n), 
in  terms  of  the  ultimate  strength,  R(0) , the  number  of  load  cycle, 
n,  and  the  maximum  stress,  S.  The  shape  of  strength  degradation 
given  in  Eq.  4 depends  on  the  value  of  c (Refs.  1,2). 

It  is  reasonable  (Refs.  1-7)  to  assume  that  the  statistical 
distribution  of  the  ultimate  strength,  R(0),  follows  a two-parameter 
Weiball  distribution. 


in  which  a is  the  shape  parameter  and  6 is  the  scale  parameter 
or  the  characteristic  strength. 

It  is  assumed  that  fatigue  failure  occurs  as  soon  as  the  resi- 
dual strength  R(n)  is  reduced  to  the  value  of  maximum  stress  level 
S.  Let  N be  the  number  of  cycles  at  which  fatigue  failure  occurs 
under  the  maximum  cyclic  stress  S.  Then,  at  the  moment  of  fatigue 
fracture,  R(n)  = S and  n = N.  The  number  of  cycles  to  fatigue 
failure  N can  be  obtained  from  Eq.  4 by  imposing  the  condition  of 
fracture,  i.e.,  R(n)  = S,  n = N,  as  follows;  = R‘^(0)  - f(S)N. 
Hence,  one  obtains 

N = [r^(0)  - /f(S)  (6) 


Since  the  ultimate  strength,  R(0),  is  a statistical  variable, 
it  follows  from  Eq.  6 that  the  number  of  cycles  to  fatigue  failure, 
N,  is  also  a statistical  variable.  The  statistical  distribution 
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of  R(0)  given  by  Eq.  5 through  the  transformation  of  Eq.  6 as 
follows ; 


F^Cn)  = P 


N < n] 

= P[j 

R(0)  < 

jnf(s) 

= P 


Substitution  of  Eq.  5 into  Eq. ' 7 yields 


F^(n)  = 1 - exp 


(7) 


r n + .S‘^/f(S\"| 

a/c  j 

L B<^/f(S)  -• 

) 

CS) 


It  is  observed' from  Eq.  8 that  the  statistical  distribution 
of  fatigue  life  N (number  of  cycles  to  fatigue  failure)  follows 
a three  parameter  Weibull  distribution  with  a lower  bound  at 
Wq  = -S^/f(S). 

For  the  sake  of  clarity,  two  different  situations  will  be 
considered  separately  in  the  following,  although  they  lead  to 
identical  residual  strength  degradation  model. 

(i)  Low  stress  amplitude  fatigue 

In  general,  when  we  are  interested  in  the  region  where  the 
maximum  stress  amplitude  S is  small  such  that  (S/6)^<<  1,  i.e., 
the  number  of  cycles  to  failure  is  high,  Eq.  8 can  be  approximated 
in  the  following; 

|a/c 

7f(sU  ) (9) 


Fjj(n)  = 1 - exp 


indicating  that  if  the  ratio  of  the  stress  level  S to  the  charac- 
teristic ultimate  strength  6 is  such  that  (S/B)^  <<  1,  the  statis- 
tical distribution  of  fatigue  life  N is  approximately  a two-parameter 
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Weibull  distribution  with  the  shape  parameter  g/c  and  the  scale 
parameter  (characteristic  life)  6^/f(S). 

Supposing  that  the  constant  amplitude  fatigue  data  is 
available,  one  obtains  a classical  S-N*  curve  by  connecting  the 
characteristic  lives  associated  with  each  stress  level  S.  The 

S-N*  curve  can  be  expressed  classically  as 

KS^N*  = 1 (10) 

in  which  K and  b are  constants  and  N*  is  the  characteristic  life 
associated  with  a maximum  stress  level  S.  Eq.  10  can  be  written 
as 

N*  = (KS'^)*^  (11) 

Comparison  of  Eqs.  9 and  11  leads  to  the  result, 
g^/f(s)  = (KS^)'^.  Hence,  the  functional  form  of  f(S)  is  obtained 
as 

f(S)  = (12) 

Substituting  Eq.  12  into  Eq.  4,  one  obtains 

R^(n)  = R^^(O)  - 6‘^KS^n  (13) 

in  which  c,  b and  K are  constants  to  be  determined  from  experi- 
mental data.  Eq.  13  is  the  equation  describing  the  residual 
strength  degradation  under  constant  amplitude  cyclic  loads. 

If  Eq.  12  is  substituted  into  Eq.  3,  one  obtains  a general 
equation 
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R^^Cnj)  = R^^Cno)  - (n^  - n^) 


(14) 


It  is  mentioned  that  the  model  presented  herein,  Eqs.  13,  and 
14,  is  similar  to  the  one  proposed  by  Halpin  and  Waddoups  (Refs. 

3-4).  However,  they  derived  the  strength  degradation  based  on 
the  assumption  of  crack  propagation.  The  model  derived  herein 
is  based  on  the  assumption  that  the  residual  strength  is  a mono- 
tonically  decreasing  function  of  load  cycles  n. 

(ii)  High  stress  amplitude  fatigue 

For  composite  laminates  where  the  endurance  limit  is  rather 
high,  fatigue  failure  will  not  occur  unless  the  maximum  stress  S 
is  very  high.  Under  this  circumstance,  (S/3)^  may  not  be  negli- 
gible, although  c is  usually  large,  e.g.,  c = 10  for  Gr/E  composites. 

Letting 

Y = (S/B)'";  Nj  = 3^/f(s)  (15) 

and  introducing  a new  random  variable  N,  indicating  a shifted 
fatigue  life  by  an  amount  of  yN^  to  the  right,  i.e., 

N = N + yNj  (16) 

one  obtains  from  Eq.  8 the  statistical  distribution  of  the  shifted 
fatigue  life  N as  follows; 

Fj;j(n)  = P [n  5 n]  = 1 - exp|  - (n/N^)“/^  j (17) 

It  follows  from  Eq.  17  that  the  statistical  distribution  of 
the  shifted  fatigue  life  N is  a two  parameter  Weibull  distribution 
with  a shape  parameter  a/c  and  a characteristic  life  of  N^. 
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Supposing  that  the  fatigue  data  is  available  for  N,  one 

I 

obtains  a classical  S-Nj^,  curve  by  connecting  the  characterstic 
' lives  associated  with  each  stress  level  S.  Such  a S-N^  curve  can 

be  approximated  by 


KS^Nj^  = 1 


(18) 


in  which  K and  b are  constants  which  are  different  from  K and 
b in  Eqs.  10-14. 


Substitution  of  appearing  in  Eq.  18  into  Eq.  15  leads 
to  the  functional  form  for  f(S), 

fCS)  = Ke'^s'^  (19) 

Substituting  Eq.  19  into  Eq.  4,  one  obtains 

R*^(n)  = r‘^(0)  - K6*^S^n  (20) 

Eq.  20  is  identical  to  Eq.  13  except  that  the  meaning  of 
the  constants  K and  b are  different.  This  indicates  that  the  form 
of  the  fatigue  residual  strength  degradation  model  of  Eq.  20  can 
be  used  for  both  the  low  stress  amplitudeandthe  high  stressamplitude  fatigue. 

It  is  very  important  to  emphasize  that  in  the  process  of 
deriving  the  residual  strength  degradation  model,  Eq.  20  or  13, 
parameters,  c,  b and  K are  assumed  to  be  determined  from  the  stress- 
characteristic  life  S-N  curve  of  statistical  fatigue  data  (see 
Eqs.  10  and  18).  The  generation  of  statistical  data  for  Eqs.  10 
and  18  is  extremely  tedious  and  should  not  be  attempted.  It  only 
serves  as  a vehicle  for  us  to  derive  the  residual  strength  degrada- 
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tion  model  only.  More  efficient  analysis  technique  and  test  pro- 
cedure will  be  presented  in  the  following,  such  that  the  number  of 
test  data  required  for  the  determination  of  parameter  values  c,  b 
and  K can  be  minimized. 

Once  the  parameter  values  K,  b and  c have  been  determined 
from  experimental  data,  Eq.  20  can  be  used  easily  to  find  the 
statistical  distribution  of  the  residual  strength,  R(n),  after 
n cycles,  and  that  of  the  fatigue  life  N for  any  value  of  maximum 
cyclic  stress  S. 

The  statistical  distribution  of  the  residual  strength  R(n) 
after  n fatigue  cycles  can  be  obtained  from  the  statistical  dis- 
tribution of  R(0)  given  by  Eq.  5 through  the  transformation  of 
Eq.  20  as  follows; 


in  which  = 1 - exp|  - +KB‘^S^n)/6^  j “'^^j  is  the 

probability  that  the  specimen  will  fracture  before  n cycles,  i.e., 
the  probability  that  the  specimen  will  not  survive  n stress  cycles, 
Note  that  there  is  a discontinuity  in  the  distribution  function, 
FR(n)(x),  at  X = 0 , and  hence  there  is  a dirac  delta  at  x = 0 
for  the  probability  density  function. 


The  number  of  cycles  N to  fatigue  failure  can  be  obtained 
from  Eq.  20  by  imposing  the  condition  of  fracture,  i.e.,  n = N, 
R(n)  = S,  as  follows; 


N = [r^(0)  - (23) 

The  statistical  distribution  of  fatigue  life  N can  be  ob- 
tained from  that  of  the  static  strength  R(0)  given  by  Eq.  5 
through  the  transformation  of  Eq.  23  as  follows: 


Ff^(n) 


P [^N  < n] 
1 - exp  I 
1 - exp  I 


= P TrCO)  < (nK6‘^S^  + 

[ fnKgCsf’  ^ S‘^)/S‘^J  I 
r n ^ (S*^/K6CS^)  1 ) 

L 1/KSb  J ' 


(24a) 


Eq.  24a  indicates  that  the  fatigue  life  distribution  is  a 
three  parameter  Weibull  distribution  with  a lower  bound  at 
KS*^.  This  comes  from  the  fact  that  some  of  the  specimens  will 
fail  within  one  cycle,  when  their  ultimate  strengths  are  smaller 
than  S.  As  a result,  it  should  be  written  as 
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I 
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F^(n)  = 0 


; n < 1 


= 1 - exp 


, a/c 

i-[ I-' 


(24b) 


where  is  the  probability  of  failure  in  one  cycle,  i.e.,  the 

probability  that  the  ultimate  strength  is  smaller  than  the  maximum 
applied  stress  S.  Again,  there  is  a discontinuity  in  the  distri- 
bution function  Fj,j(n)  at  n = 1,  and  hence  a dirac  delta  at  n = 1 
in  the  probability  density  function. 

(2)  Test  Plan  and  Analysis  Technique  for  the  Determination  of 
Parameter  Values. 


Economical  test  plan  and  efficient  analysis  technique  are 
presented  herein  to  minimize  the  number  of  test  data  required  for 
the  determination  of  parameter  values  c,  b and  K appearing  in 
the  residual  strength  degradation  model.  Let  Xj,  X2»  ..., 
be  static  ultimate  strengths  of  a set  of  m specimens,  referred 
to  as  static  specimens.  The  first  three  central  moments  of  the 
static  strength  data  are  given  as  follows; 

m 

•"I  = iii  ^i 


m 

Z 

i=l 


m 


3 


1 

m 


m 

Z 

i = l 


(x^-mj)' 


(25) 


Note  that  mj^  is  the 
deviation  and  m|^^/  m^ 


sample  mean,  /in^  is  the  sample  standard 
is  the  sample  skewness 


A set  of  .J  specimens  are  subjected  to  constant  amplitude 
fatigue  at  various  maximum  stress  level  S,  referred  to  as  fatigue 
spec imens . Some  specimens  will  fracture  under  fatigue  within 
10^  or  10^  cycles,  and  some  don't.  When  the  specimen  does  not 
fail,  the  fatigue  test  is  stopped  and  its  residual  strength  is 
measured.  Furthermore,  the  residual  strength  of  the  specimen 
at  the  instant  of  fracture  in  fatigue  is  equal  to  the  maximum 
cyclic  stress  S.  As  a result,  a set  of  J data  points  of  resi- 
dual strengths,  denoted  by  (Rj^ , , n^j  , (R2»  S2,  n2) 

(Rj,  Sj,  nj) , is  obtained  from  fatigue  specimens.  Note  that  the  ith 
data  point  (Rj^,  n^^)  , indicates  that  the  ith  specimen 

is  subjected  to  a maximum  cyclic  stress  for  n^  cycles  and 
its  residual  strength  is  R^^.  If  the  specimen  fracture  at  nj^  cycles, 
then  Rj^  = Sj^. 

Now  the  residual  strength  R^(i  = 1,2 J)  of  the  fatigue 

specimens  can  be  converted  into  (or  recovered  to)  the  equivalent 
ultimate  strength,  denoted  by  R^(0),  through  the  strength  degrada- 
tion model  of  Eq.  20,  i.e., 

R.(0)  = [r?  + KB‘^S^n-1  ; i = 1,2,...,  J. 

^ L ^ ^ (26) 

in  which  R^(0)  would  have  been  the  static  ultimate  strength  if 

the  ith  fatigue  specimen  were  not  subjected  to  fatigue. 

Using  Eq.  26,  one  obtains  from  fatigue  data  a set  of  J equiva- 
lent ultimated  static  strength  data,  i.e.,  Rj^(O),  R2(0)...,  Rj(0). 

The  first  three  central  moments,  denoted  by  Uj,  M2  i‘3 
follows ; 
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Ml  = 4 Z R.(0) 

i = l 1 

^2  " J j^Ri(O) 

.3-1  [«i"» 


Theoretica] ly , if  the  number  m of  static  ultimate  strength 
data  and  the  number  J of  fatigue  data  approach  infinity,  the 

statistical  distribution  of  x^(i=l,2 m)  should  be  identical 

to  the  statistical  distribution  of  R-(0)  (i=l  ,2 , . . . ,J)  , and  hence 
mj^  = Mj,  ni2  = ^2  and  = M3,  if  the  residual  strength  degra- 
dation model  (Eq.  20)  is  correct.  In  reality,  however,  since 
the  number  of  data  points  is  finite,  sampling  fluctuation  can 
not  be  avoided, and  hence  it  may  not  be  possible  to  match  the 
first  three  central  moments  to  determine  c,  b and  K.  Consequent- 
ly, c,  b and  K are  obtained  by  minimizing  the  mean  square  differ- 
ence A between  the  first  three  central  moments  i.e., 

2 2 

A = (mj^  - M^)^  + ^ ) + g3(^/jii7  ' ) 

(28) 

where  g^  and  g^  are  assigned  weighting  positive  values  to  indicate 
the  relative  importance  of  matching  the  mean,  the  standard  devia- 
tion and  the  skewness. 

It  becomes  clear  now  that  in  order  to  determine  the  parameter 
values  of  the  residual  strength  degradation  model  (Eq.  20),  i.e., 
K,  b and  c,  one  needs  to  test  a set  of  m static  data  and  a set 
of  J fatigue  data.  Since  the  statistical  dispersion  of  static 
strength  data  is  small,  m is  not  necessarily  large.  The  total 
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number  of  fatigue  data  J depends  on  the  ranges  of  fatigue  life 
and  stress  levels  which  are  of  interest  to  us.  Experience  indi- 
cates that  J ~ 30  is  sufficient  for  the  stress  range  S ^ 60%  of 
the  mean  ultimate  strength. 


Experimental  Verification 


A 2'  X 2'  panel  of  graphite/epoxy  laminate  was  curve  using 
the  Gr/E  tape  purchased  from  Hercules  AS3501.  This  panel,  referred 
to  as  panel  3,  was  cut  into  approximately  70  specimens.  The  layup 
of  the  laminates  is  (0,90,  ± 45)^  and  the  size  of  the  specimen 
is  shown  in  Fig.  1. 

Twelve  (12)  specimens  are  tested  statically  and  their  static 
ultimate  strengths  are  shown  in  Table  I.  The  two-parameter 
Weibull  distribution  (Eq.  5)  is  used  to  fit  the  data  as  shown 
by  the  solid  curve  of  Fig.  2(a).  A shape  parameter  a = 14.778 
and  a characteristic  ultimate  strength  8 = 78.109  ksi  have  been 
obtained. 


Twenty  one  (21)  specimens  are  subjected  to  fatigue  cycles 
at  various  maximum  stress  levels  S and  the  results  are  tabulated 
in  Table  II.  One  specimen.  No.  2-12,  did  not  fail  at  one  million 
cycles  and  its  residual  strength  was  measured.  In  addition,  4 
specimens  are  fatigued  at  a maximum  stress  level  72.5%  of  the 
mean  ultimate  strength  y = 75.391  ksi  up  to  30,000  cycles.  The 
fatigue  test  is  then  stopped  and  their  residual  strengths  are 
measured.  Results  are  shown  at  the  bottom  of  Table  II.  All 
the  fatigue  tests  are  performed  at  a frequency  of  20Hz  with  a 
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stress  ratio  0.1  on  an  Instron  cyclers.  Both  the  static  data 
and  the  fatigue  data  (Table  I and  II)  are  plotted  in  Fig.  3. 

A total  of  12  static  data  points  in  Table  I and  25  fatigue 
data  points  in  Table  II  are  used  to  determine  the  parameter  values 
c,  b and  K using  the  method  of  analysis  described  in  Eq.  25-28; 
with  the  result, 

C = 10.818;  b = 17.78;  K = 1.8285  x 10'^^  (29) 

As  a first  indication  of  the  goodness  of  the  residual  strength 
degradation  model  (Eq.  20),  a set  of  25  fatigue  data  (Table  II) 
is  converted  into  a set  of  equivalent  static  ultimate  strength 
data  R^(0)  (i  = 1,2,..., 25)  through  Eq.  26  using  c,  b and  K 
values  given  by  Eq.  29.  Then,  this  set  of  25  equivalent  ultimate 
strength  data  is  normalized  with  respect  to  B and  is  plotted  in 
Fig.  2(b)  as  circles  where  the  solid  curve  is  the  Weibull  distri- 
bution of  12  static  strength  data  (Table  I). 

It  can  be  observed  from  Fig.  2(b)  that  the  correlation  between 
the  fatigue  converted  ultimate  strength  data  (circles)  and  the 
static  ultimate  strength  distribution  is  outstanding.  This  is 
a necessary  condition  for  the  residual  strength  degradation  model, 
Eq.  20,  to  be  reasonable. 

Since  R(0)  and  R(n)  are  random  variables,  Eq.  20  can  also 
be  written  as 

R^(n)  = Rp(0)  - KB^S^n  (30) 
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in  which  Rp(n)  and  Rp(0)  represent,  respectively,  the  p reliability 
point.  A plot  of  Eq.  30  is  given  by  Fig,  4 along  with  25  fatigue 
data  (Table  II). 

With  the  parameter  values  of  c,  b and  K as  well  as  Eq.  20, 
it  is  possible  to  predict  the  residual  strength  distribution 
under  any  maximum  stress  level  S at  any  number  of  cycles  n by 
use  of  Eq.  22. 

To  verify  the  validity  of  the  residual  strength  degradation 
model  (Eq.  20) , a set  of  8 specimens  are  fatigue  tested  at  the 
maximum  stress  level  of  S = 0.7p  = 52.77ksi  up  to  26,000  cycles. 

One  specimen  fails  at  17,540  cycles.  The  residual  strengths  of 
the  remaining  7 surviving  specimens  are  given  in  Table  III. 

According  to  the  theoretical  model,  Eq.  22,  the  probability 
of  failure  before  26,000  cycles  is  Fp  ^^^^(S)  - 0. 1135  (see  Eq. 

22).  The  experimental  data  shows  one  failure  out  of  8 specimens, 
i.e.,  1/(8  + 1)  = 0.11,  indicating  that  the  theoretical  prediction 
is  very  good.  The  theoretical  prediction  (Eq.  22)  is  plotted  in 
Fig.  5(a)  as  a solid  curve  and  the  experimental  data  as  circles. 
Considering  the  sampling  fluctuation  due  to  small  sample  size, 
the  correlation  between  test  data  and  theoretical  prediction 
appears  to  be  reasonable. 

Furthermore  the  conditional  distribution  function  of 
the  residual  strength  R* (n)  for  surviving  specimens  can  also 
be  predicted  theoretically  using  Eq.  22  as  follows; 
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= 1 - P j^R(n)  > X I 
= 1 - P ^R(n)  > X 
= 1 

(30) 

in  which  (S)  is  given  by  Eq . 22. 

The  theoretical  prediction  (Eq.  30)  is  plotted  in  Fig.  5(b) 
as  solid  curve  and  the  experimental  data  as  circles.  Again,  the 
correlation  is  reasonable  in  view  of  small  sample  size. 

SECTION  III 
PERIODIC  PROOF  TESTS 
(1)  Experimental  Data 

A set  of  14  specimens  is  subjected  to  a proof  load  r^  = 0 . 9y  = 
67.816  ksi  prior  to  fatigue.  After  the  specimen  passes  the  initial 
proof  load  rg,  it  is  subjected  to  26,000  fatigue  cycles  at  the 
maximum  stress  level  S = 0.7p  = 52.78  ksi.  A second  proof  load 
rQ  is  further  applied  to  the  specimen  and  then  followed  by  another 
26,000  fatigue  cycles.  The  residual  strengths  of  surviving  speci- 
mens are  measured  at  this  point.  A loading  history  is  schematically 
shown  in  Fig.  6.  Results  of  residual  strengths  are  shown  in  Table 
IV. 

Note  that  there  is  a possibility  that  the  specimen  may  fail 
under  the  initial  proof  load  or  the  second  proof  load.  Test 
results  show  that  3 specimens  fail  under  the  2nd  proof  test  and 
their  residual  strengths  are  indicated  by  a star  in  Table  IV.  No 
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specimen  fails  under  the  initial  proof  test  or  by  fatigue. 

A second  set  of  11  specimens  is  subjected  to  the  same 
loading  history  as  shown  in  Fig.  6,  except  that  the  proof  load 
level  is  851  of  the  mean  static  ultimate  strength,  i.e.,  rQ  = 

0.85p  = 64.1  ksi.  Residual  strengths  of  8 surv.^ing  specimens 
are  given  in  Table  V.  Three  specimens  fail  under  the  second 
proof  test  and  their  residual  strengths  are  indicated  in  Table 
V by  a star.  No  specimen  fails  under  the  initial  proof  test 
or  by  fatigue. 

(2j  Theoretical  and  Experimental  Correlation 

Theoretical  derivation  for  the  statistical  distribution  of 
the  residual  strength  of  surviving  specimens  under  periodic  proof 
tests  and  service  load  has  been  presented  in  Ref.  1.  The  results 
for  the  particular  case  of  constant  amplitude  fatigue  is  presented 
in  the  Appendix.  Theoretical  predictions  for  testing  conditions 
associated  with  two  sets  of  experimental  data  are  presented  in 
Fig.  7 and  8 as  solid  curves,  along  with  the  experimental  data 
indicated  by  circles. 

It  can  be  observed  from  Figs.  7 and  8 that  the  correlation 
between  the  theory  of  periodic  proof  tests  and  the  experimental 
results  is  surprisingly  good  for  both  cases. 

According  to  the  theoretical  predictions  derived  in  the 
appendix,  the  following  conclusions  are  made  for  the  first  set 
of  experimental  data  (Table  IV). 
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(i)  If  the  specimen  passes  a proof  load  = 0.9p  = 67.852  ksi, 
it  will  survive  at  least  26,750  fatigue  cycles  at  the  maximum  stress 
level  S = 0.7p  = 52.7  ksi.  This  can  by  computed  easily  from  Eq. 

23  by  setting  R(0)  = r^.  This  conclusion  is  confirmed  by  the 
experimental  data  that  none  of  the  14  specimens  fails  in  fatigue,  i.e., 
when  the  specimen  passes  either  the  initial  proof  test  or  the  second 
proof  test,  it  does  not  fail  in  the  subsequent  26,000  cycles  of 
fatigue  load.  This  is  an  essential  and  important  demonstration 
that  the  periodic  proof  tests  not  only  sensors  unnotched  specimens 
that  are  weak  in  static  strength,  but  also  screens  out  specimens 
having  weak  fatigue  strength,  thus  guarranteeing  a minimum  fatigue 
life  for  specimens  when  they  pass  certain  level  of  proof  load. 

This  is  an  essential  advantage  for  the  application  of  periodic 
proof  tests  as  advocated  in  Refs.  1-4. 

(ii)  The  probability  of  failure  of  a specimen  under  the  ini- 
tial proof  load  rg  = 0.9y  is  computed  from  Eq.  A - 4 in  the  Appen- 
dix as  0.117.  The  probability  of  failure  of  a specimen  under 

the  second  proof  test  is  computed  from  Eq.  A - 7 in  the  Appendix 
as  0.161.  As  a result,  the  total  probability  that  a specimen  will 
fail  under  two  proof  tests  is  0.278.  Consequently,  for  the  14 
specimens  tested,  an  average  of  14  x 0.278  = 3.9  specimens  will 
be  destroyed  by  proof  tests.  The  experimental  results  show  3 
specimens  failed  under  the  second  proof  test.  Considering  the 
sampling  fluctuating  associated  with  small  sample  size,  the  corre- 
lation between  the  theory  and  the  experimental  data  is  reasonable. 
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(iii)  The  statistical  distribution  of  the  residual  strength 
of  the  specimens  surviving  2 proof  tests  and  two  service  intervals, 
where  each  interval  is  26,000  fatigue  cycles  at  S = 0.7y,  has 
been  derived  in  Eq.  A-6  of  the  Appendix.  This  theoretical  distri- 
bution is  plotted  in  Fig.  7 as  solid  curve.  Also  plotted  in  Fig. 

7 as  circles  are  the  residual  strengths  of  the  experimental  data 
(Table  V).  It  can  be  observed  that  the  correlation  between  the 
theoretical  prediction  and  experimental  results  is  excellent. 

(iiii)  The  fact  that  good  correlation  exists  between  the 
theory  of  periodic  proof  tests  and  experimental  data  indicates 
that  both  the  residual  strength  degradation  model  and  the  theory 
of  periodic  proof  tests  are  reasonable  and  can  be  used  with  cer- 
tain degree  of  confidence.  No  good  correlation  can  be  obtained, 
if  either  the  residual  strength  degradation  model  or  the  theory 
of  periodic  proof  tests  is  unreasonable. 

The  conclusion  for  the  second  set  of  experimental  data  are 
given  in  the  following: 

(i)  According  to  Eq.  23,  the  specimen  will  sur»/^ive  at  least 
14,000  fatigue  cycles  at  S = 0.7u  = 52.78  ksi  after  passing  a 
proof  load  r^  = 0.85y  = 64.1  ksi.  Since  the  specimens  are  sub- 
jected to  two  service  intervals  each  with  26,000  fatigue  cycles, 
there  is  a slight  probability  that  the  specimen  may  fail  by 
fatigue.  The  probability  of  failure  under  initial  proof  test 
is  computed  from  Eq.  A - 4 as  0.0525.  The  probabilities  of  failure 
(a)  in  the  first  service  interval,  (b)  under  the  second  proof 
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test  and  (c)  in  the  second  service  interval,  are  computed, 
respectively,  as  = 0.064,  = 0.078  and  ?2  ~ 0.073  (see 

Eqs.  A-11  to  A-13  of  the  Appendix).  As  a result,  the  total 
probability  that  the  specimen  will  not  survive  up  to  the  end  of 
the  second  service  interval  is  0.2675.  Hence,  the  average 
number  of  specimens  expected  to  fail  either  in  fatigue  or 
under  proof  tests  is  0.2675  x 11  = 2.94  specimens.  The  exper- 
imental results  show  (Table  III)  that  3 specimens  fail  before 
the  end  of  the  second  service  interval.  The  correlation  between 
the  theory  and  experimental  results  is  good. 

(ii)  The  theoretical  prediction  for  the  residual  strength 
distribution  of  surviving  specimens  has  been  derived  in  Eq. 

A - 17  of  the  Appendix,  and  is  plotted  as  solid  curve  in  Fig.  8. 
Also  plotted  in  Fig.  8 as  circles  are  the  results  of  experimental 
data.  Again,  the  correlation  between  the  theory  of  periodic 
proof  tests  and  experimental  data  is  excellent. 

(3)  Damage  Due  to  Proof  Tests 

It  has  been  a concern  that  proof  tests  may  damage  composites 
in  terms  of  the  reduction  of  both  the  static  strength  and  the 
fatigue  strength  as  discussed  in  Refs. l-4.Damage  in  static  strength 
may  come  from  the  fact  that  under  one  cycle  of  large  proof  load, 
the  existing  crack  may  increase  its  size  thus  reducing  the  static 
strength.  It  has  been  shown  experimentally  that  such  a damage 
is  negligible  for  Gr/E  tt/4  laminates  (Ref.  9). 

Damage  in  fatigue  strength  may  come  from  the  initiation  of 
cracks,  e.g.,  first  ply  failure,  due  to  a high  proof  load.  Our 
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theoretical  predictions  derived  in  the  Appendix  are  based  on  the 
assumption  that  no  damage  in  fatigue  strength  is  produced  by 
proof  tests.  The  fact  that  excellent  correlation  exists  between 
experimental  results  and  theoretical  predictions  indicates 
that  damage  is  negligible  for  the  proof  load  r^  up  to  90%  of 
the  mean  ultimate  strength. 

In  fact,  the  ideal  experimental  testing  condition  for  the 
verification  of  damage  in  fatigue  strength  is  by  use  of  cyclic 
loading  that  is  below  the  first  ply  failure  stress.  It  should 
be  noticed,  however,  that  the  endurance  limit  for  Gr/E-n/4  laminates 
is  above  the  first  ply  failure  stress.  As  a result,  cyclic 
loading  below  the  first  ply  failure  stress  is  of  no  concern  in 
design . 

SECTION  IV 

Conclusion  and  Recommendation 

A theoretical  fatigue  residual  strength  degradation  model 
has  been  derived  based  on  the  assumption  that  the  residual 
strength  decreases  monotonically . Such  a model  is  applicable 
to  unnotched  laminates  as  demonstrated  herein,  as  well  as  adhe- 
sive or  mechanically  fastened  joints  (Refs.  6-8).  It  is  expected 
to  be  applicable  to  laminates  with  matrix  dominated  failure  mode. 
This  type  of  residual  strength  degradation  model  can  also  be  ex- 
tended to  the  case  where  the  loads  are  random  in  nature,  re- 
sembling the  real  service  loads. 
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For  notched  laminates,  it  has  been  observed  that  the  resi- 
dual strength  increases  initially  and  then  decreases  eventually. 
Under  such  a situation,  it  is  recommended  that  different  theore- 
tical model  should  be  developed.  The  model  should  be  such  that 
it  approaches  the  strength  degradation  model  described  herein 
as  the  notched  size  approaches  zero  in  a limiting  case. 


We  have  demonstrated  experimentally  the  validity  of  both 
the  derived  residual  strength  degradation  model  and  the  theory 
of  periodic  proof  tests  using  a rather  small  sample  size.  The 
reason  for  using  such  a small  sample  size  is  due  to  the  limita- 
tion that  we  have  only  one  2'  x 2'  panel  with  square  end  taps 
(Fig.  1)  available.  The  panel  is  cut  into  70  specimens.  According 
to  our  experience,  end  tap  failure  can  be  avoided  by  using  equare 
end  taps.  All  the  data  reported  herein  are  not  tap  failure.  It 
is  recommended  that  large  sample  size  should  be  used  later  in 
order  to  reduce  sampling  fluctuation  due  to  the  statistical 
nature  of  the  fatigue  behavior  of  composite  laminates. 


We  have  specimens  from  one  panel,  referred  to  as  panel  1, 
with  ramped  end  tap.  Unfortunately,  the  specimen  with  this  type 
of  end  tap  results  in  end  tap  failure  under  fatigue  loads.  Limited 
amount  of  fatigue  data  (end  tap  failure)  and  static  ultimate 
strength  data  have  also  been  generated  and  they  also  appear  to 
fit  the  residual  strength  degradation  model  described  herein,  but 
with  quite  different  parameter  values  c,  b and  K.  In  other  words, 
the  shape  of  the  residual  strength  degradation  is  different  from 
data  without  end  tap  failure. 
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It  has  been  observed  that  there  may  be  aging  effect  for 
the  strength  of  composite  laminates  due  to  residual  stress 
relaxation.  Sould  this  phenomenon  occur,  the  composite  strength 
would  increase  with  time.  It  should  be  mentioned  that  all  the 
data  presented  in  this  paper  were  generated  within  a one  and 
a half  month  period.  It  is,  therefore,  reasonable  to  expect 
that  our  data  should  have  minimal  aging  effect. 


TABLE  1.  ULTIMATE  STRENGTH  Gr/E  [ 0 , 90 , ±45]^,  PANEL  3 


Specimen  No. 

Ultimate  Strength 
ksi 

1-3 

63.152 

1 

2 

1-2 

66.312 

3 

5-2 

71.900 

4 

1-1 

11. Ill 

5 

4-1 

11. tit 

6 

5-1 

75.050 

7 

3-2 

77.743 

8 

2-3 

78.316 

1 

! 9 

2-2 

80.052 

10 

4-2 

81.324 

11 

2-1 

81.742 

12 

3-1 

84.154 

Average 

75.39  i 

Standard  Deviation 

6.12 

Coefficient  of  Variation  0.08 

a 

14.78 

6 

78.1 
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TABLE  II.  FATIGUE  AND  RESIDUAL  STRENGTH  DATA  Gr/E[0,90,  ±4S 

PANEL  3 


Maximum  Stress 
Level  S in 
1 of  u 

No.  of  Cycles 
to  Fracture 

85 

1,650 

85 

1,950 

85 

1,320 

80 

2,050 

75 

50,980 

75 

6,480 

70 

155,000 

70 

228,500 

70 

88,000 

70 

117,580 

70 

228,700 

70 

221,200 

70 

310,000 

70 

18,790 

70 

3,840 

67 

161,000 

67 

110,000 

65 

523,500 

65 

863,200 

62 

1,346,300 

65 

1,007,000* 

72.5 

30,000* 

72.5 

30,000* 

72.5 

30,000* 

72.5 

30,000* 

Strength 


64.012 

64.012 

64.012 

60.246 

56.481 

56.481 

52.716 

52.716 

52.716 

52.716 

52.716 

52.716 

52.716 

52.716 

52.716 
50.456 
50.456 


48.950 

46.691 

60.201 

61.455 

64.544 


75.39  ksi  = mean  ultimate  strength. 

Specimen  does  not  fracture  in  fatigue. 

The  residual  strength  at  fatigue  fracture  is  equal  to  S 


TABLE  III.  RESIDUAL  STRENGTH  AFTER  26,000  FATIGUE 
CYCLES  AT  S = 0.7u 


Specimen  No. 


Residual  Strength 
(ksi) 


2-10 

3- 9 

4- 7 


8 

9 

10 

11 

12 

13 


1-10 

2-11 

3- 12 

4- 10 

1-16 

2-13 

2- 14 
1-13 
1-17 

3- 3 


5-14 


79.68 

63.61 

80.17 

77.24 
75.44 
82.03 
80.98 

71.24 
76.31 
67.28 
69.71 
66.10* 
62.10* 
63.59* 
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Specimen  fails  under  the  second  proof  test 


TABLE  V.  RESIDUAL  STRENGTH  AFTER  TWO  SERVICE  INTERVALS; 


Tq  = 0 . 85u  = 64.1  ksi 


Specimen  No. 

Residual  Strength 
(ksi) 

1 

1-4 

58.34 

2 

4-11 

67.20 

3 

2-14 

69.74 

4 

1-6 

73.15 

5 

2-5 

75.43 

6 

3-6 

76.02 

7 

2-6 

78.00 

8 

3-11 

80.05 

9 

5-3 

57.73* 

10 

5-5 

62.74* 

11 

4-5 

63.80* 

* Specimen  fails  under  the  second  proof  test 


DISTRIBUTION  FUNCTION 
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STATIC  ULTIMATE  STRENGTH  IN  p CONVERTED  ULTIMATE  STRENGTH 


OF  CYCLES 


NO.  OF  CYCLES 


RESIDUAL  STRENGTH  IN/A 


■f. 


f' 

I 


FIG.  4.(d)  residual  STRENGTH  DEGRADATION  UNDER  MAXIMUM  FATIGUE  STRESS 
LEVEL  S = 0.725P 

36 


NO.  OF  CYCLES 


NUMBER  OF  CYCLES 


90% 


NUMBER  OF  CYCLES 


DISTRIBUTION  FUNCTION 


I 

'P 


50  S 60  70  80  90 


RESIDUAL  STRENGTH  AFTER  26,000 
FATIGUE  CYCLE,  KSI 

FIG.  5.:  STATISTICAL  DISTRIBUTION  OF  RESIDUAL  STRENGTH; 

(a)  unconditional  distribution 

41 


CONDITIONAL  DISTRIBUTION  FUNCTION 


THEORETICAL 


RESIDUAL  STRENGTH  AFTER  26,000 
FATIGUE  CYCLES  FOR  SURVIVING 
SPECIMENS,  KSI 

FIG.  5.:  STATISTICAL  DISTRIBUTION  OF  RESIDUAL  STRENGTH: 

(b)  CONDITIONAL  DISTRIBUTION  FOR  SURVIVING  SPECIMENS 


STRESS  LEVEL  IN  MEAN  STRENGTH  fl 


DISTRIBUTION  FUNCTION 


ll 


60  70  80 

RESIDUAL  STRENGTH  R(2T  ) IN  ksi 


FIG.  7,;  STATISTICAL  DISTRIBUTION  OF  RESIDUAL  STRENGTH  UNDER 
PERIODIC  PROOF  TESTS  AND  FATIGUE  LOADS 
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DISTRIBUTION  FUNCTION  ^ R (2T-I*^* 


FIG. 


RESIDUAL  STRENGTH  R(2T  ) IN  ksi 

8.:  STATISTICAL  DISTRIBUTION  OF  RESIDUAL  STRENGTH  UNDER 
PERIODIC  PROOF  TESTS  AND  FATIGUE  LOADS 
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APPENDIX 


PROBABILITY  OF  FAILURE  UNDER  PROOF  TESTS 


Based  on  Eq.  1,  the  residual  strength  degradation  under 
ramdom  service  loads  has  been  derived  in  Ref.  1 (Eq.  4)  as 


R'^Cnj^)  = R^^CHq)  - 4>{n^  - n^) 


(A-1) 


For  n^  = n,  ng  = 0,  Eq.  A-1  becomes 


R^^Cn)  = r'^(O)  - <fn 


(A-2) 


For  the  particular  case  of  constant  amplitude  fatigue,  it 
can  be  observed  from  Eq.  20  and  Eq.  A-2  that 


<D  = b'^KS^ 


(A-3) 


Under  the  initial  proof  test  prior  to  service,  the  probabil- 
ity of  failure,  denoted  by  B„ , is  obtained  from  Eq.  5 as 


Bg  = P R(0)  < rg 


1 - exp  { - (rf^/B) 


(A-4) 


After  the  initial  proof  test,  the  statistical  distribution 
of  surviving  specimens  is  truncated  at  rg  and  can  be  derived  as  * 
(see  Ref.  1 , Eq.  2)  , 

^R(0+)^^^  = 1 - exp  (rg/B)“  - (x/B)“  ; x > rg 

(A-5) 


*In  Ref.  1,  the  ultimate  strength  of  virgin  specimen  is  denoted 
by  Rq,  while  R(0)  represents  the  ultimate  strength  of  specimens 
passing  the  initial  proof  test.  In  the  present  paper,  R(0)  and 
R(0+)  are,  respectively,  identical  to  Rg  and  R(0)  in  Ref.  1. 


With  the  assumption  of  negligible  failure  probability  in 
service,  the  statistical  distribution  of  the  residual  strength 
after  the  second  service  interval  has  been  derived  as  (see  Ref. 
1,  Eq.  1-7), 

a/c 
(A-6) 

in  which  T is  the  number  of  cycles  in  one  service  interval 
(26,000  cycles)  and  $ is  given  by  Eq.  A-3. 

The  probability  of  failure  under  the  second  proof  test  is 
given  by  (Ref.  1,  Eq.  1-5), 

1 - B*  = 1 - exp  [(rp/B)"  - [(rp  + 

Equations  given  above  are  applicable  to  the  case  where 
rp  = 0.9y,  i.e.,  data  in  Table  IV,  since  the  probability  of 
failure  in  each  service  interval  is  zero.  For  the  case  where 
rp  = 0.85y, the  probability  of  failure  in  each  service  interval 
is  not  negligible,  because  the  specimen  passing  the  proof  test 
can  be  guarranteed  to  survive  at  least  14,000  cycles  only. 
Consequently,  the  theoretical  predictions  for  the  data  presented 
in  Table  V have  to  be  rederived  as  follows; 

Let  r^^,  and  r^  be  the  static  ultimate  strength  above 
which  the  specimen  can  survive,  respectively,  the  first  service 
interval  (26,000  cycles),  the  second  proof  test,  and  second 
service  interval. 


^^R(2T-) 


1 - exp 


I’D  + <I>T 


.a/c 


oC 

p 


x*~  + 2$T 
BC 


Vj  is  computed  by  setting  R(0)  = and  R(n)  = S = 0.7u  in 


equation  20  for  n = 26,000 


1/c 


= 67.684  ksi 


fA-8) 


r.,  is  computed  by  setting  R(0)  = r.,  and  R(n)  = = 0.85m 

in  F.q.  20  for  n = 26,000  and  S = 0.7m 


1/c 


= 70.21  ksi 


fA-9) 


Finally,  is  computed  by  setting  R(0)  = r^,  and  R(nj^)  = 
S = 0.7m  in  Fq.  20  for  n^  = 52,000 


r j ^ + B'^KS^n^  j 


1/c 


71.93  ksi 


(A-10) 


The  probability  of  failure  in  the  first  service  interval, 
denoted  by  , is  the  percentage  of  specimens  having  ultimate 
strength  between  r^  and  r^,  i.e., 

( 

. a 


b ■ - 1 - «>'P 


(r^/8)“  - (r^/6)' 


= 0.064 
(A-11) 


The  probability  of  failure  under  the  second  proof  test, 
denoted  by  , is  the  percentage  of  specimens  having  ultimate 
strength  between  r^  and  r^ , i.e., 


4 ^R(0+)*^’’2^  ' ^R(0  + )^’^1^  ~ 


(A-12) 


in  which  is  given  by  Eq.  A-5. 

Similarily,  the  probability  of  failure  in  the  second  service 
interval,  denoted  by  ?2  is 


-48- 


#*V  tM'  « A *(k‘  »•  - 


V*v  • 'v  ' 


(A-13) 


r 


t 

I « 


( !• 


i 

‘ ■V’ 


I „ J 
I 

I 

t. 


^2  ^Ri0  + )^^3^  ■ ^Rro  + )^^2^  - (A-13) 

As  a result,  all  the  specimens  with  ultimate  strength  less 
than  r^  will  not  survive  up  to  the  end  of  the  second  service 
interval.  The  statistical  distribution  of  the  ultimate  strength, 
R(0),  of  the  surviving  specimens  is  therefore  truncated  at  r^. 


^R(C)  ^ 

= 1 - P 


R(0+)  > X I R(0+)  > r 


R (0 


+ ) >xj  / P R(0  + ) > r3  j; 


X > r3 

(A-14) 


Substitution  of  Eq.  A-5  into  Eq.  A-14  leads  to  the  following 

\ 

^RCO)*^"^^  = 1 - exp  j (r3/3)°‘-  Cx/S)“[-,x  > r3  (A-15) 

The  ultimate  strength,  R(0) , of  the  surviving  specimens  is  sub- 
jected to  degradation  following  Eq.  20,  i.e.. 


R^(n)  = R'^CO)  - B^KS^n 


(A-16) 


in  which  n = 52,000  and  R(n)  is  the  residual  strength  at  the 
end  of  the  second  service  interval.  The  statistical  distribution 
for  R(n)  is  obtained  from  that  of  R(0)  given  by  Eq.  A-15  through 
the  transformation  of  Eq.  A-16  as  follows; 

a/c; 

BCKS% 

( 

(A-17) 


^R(n) (^)  - 1 - 


in  which  r^  = 71,93  ksi  and  n = 52,000. 
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